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ABSTRACT 

In  this  paper  we  further  explore  a  class  of  high  order  TVD  (total  variation  diminishing) 
Runge-Kutta  time  discretization  initialized  in  [12],  suitable  for  solving  hyperbolic  conser¬ 
vation  laws  with  stable  spatial  discretizations.  We  illustrate  with  numerical  examples  that 
non-TVD  but  linearly  stable  Runge-Kutta  time  discretization  can  generate  oscillations  even 
for  TVD  (total  variation  diminishing)  spatial  discretization,  verifying  the  claim  that  TVD 
Runge-Kutta  methods  are  important  for  such  applications.  We  then  explore  the  issue  of 
optimal  TVD  Runge-Kutta  methods  for  second,  third  and  fourth  order,  and  for  low  storage 
Runge-Kutta  methods. 
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1  Introduction 


In  this  paper  we  further  explore  a  class  of  high  order  TVD  (total  variation  diminishing) 
Runge-Kutta  time  discretization  initialized  in  [12].  For  related  work  of  multi-step  type  see 
[11].  The  method  is  used  to  solve  a  system  of  ODEs: 

ut  =  L{u)  (1.1) 

with  suitable  initial  conditions,  resulting  from  a  method  of  lines  approximation  to  a  hyper¬ 
bolic  conservation  law: 

Ut  =  -f{u)x  (1-2) 

where  the  spatial  derivative  f{u)x  is  approximated  by  a  TVD  finite  difference  or  finite  element 
approximation  (e.g.  [4],  [8],  [13],  [2]),  denoted  by  -L{u),  which  has  the  property  that  the 
total  variation  of  the  numerical  solution: 


does  not  increase 


TV{u)  =  Y,  |wj+i  - 

3 


rv(u"+^)  <  ry(u") 


for  a  first  order  in  time  Euler  forward  stepping: 


(1.3) 


(1.4) 


=  ti"  +  AtX(u”)  (1.5) 

under  suitable  restriction  on  At: 

At  <  Ati.  (1.6) 

The  objective  of  the  high  order  TVD  Runge-Kutta  time  discretization,  is  to  maintain  the 
TVD  property  (1.4)  while  achieving  higher  order  accuracy  in  time,  perhaps  with  a  different 
time  step  restriction  than  (1.6): 

At  <  cAti.  (1.7) 

where  c  is  termed  CFL  coefficient  for  the  high  order  time  discretization. 
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The  TVD  high  order  time  discretization  is  useful  not  only  for  TVD  spatial  discretizations, 
but  also  for  TVB  (total  variation  bounded)  (e.g.  [10]),  or  ENO  (Essentially  Non-Oscillatory) 
(e.g.  [5],  [12]),  or  other  types  of  spatial  discretizations  for  hyperbolic  problems.  It  maintains 
stability  in  whatever  norm,  of  the  Euler  forward  first  order  time  stepping,  for  the  high  order 
time  discretization,  under  the  time  step  restriction  (1.7).  For  example,  if  it  is  used  for  multi 
space  dimensional  scalar  conservation  laws,  for  which  TVD  is  not  possible  but  maximum 
norm  stability  can  be  maintained  for  high  order  spatial  discretizations  plus  Euler  forward 
time  stepping  (e.g.  [3]),  then  the  same  maximum  norm  stability  can  be  maintained  if  TVD 
high  order  time  discretization  is  used.  As  another  example,  if  an  entropy  inequality  can  be 
proved  for  the  Euler  forward,  then  the  same  entropy  inequality  is  valid  under  a  high  order 
TVD  time  discretization. 

In  [12],  a  general  Runge-Kutta  method  for  (1.1)  is  written  in  the  form: 

+  i  =  (1.8) 

k=:0 

Clearly,  if  all  the  coefficients  are  nonnegative  aik  >  0,  ^ik  >  0,  then  (1.8)  is  just  a  convex 
combination  of  Euler  forward  operators,  with  At  replaced  by  since  by  consistency 

E1=o  =  1-  We  thus  have 

Lemma  1.1.  [12]  The  Runge-Kutta  method  (1.8)  is  TVD  under  the  CFL  coefficient  (1.7): 

c  =  min^,  (1-9) 

/3ik 

provided  that  aik  >  0,  ^ik  >  0. 

□ 

In  [12],  schemes  up  to  third  order  were  found  to  satisfy  the  conditions  in  Lemma  1.1  with 
CFL  coefficient  equal  to  1. 

If  we  only  have  aik  >  0  where  I3ik  might  be  negative,  we  need  to  introduce  an  adjoint 
operator  L.  The  requirement  for  L  is  that  it  approximates  the  same  spatial  derivative(s)  as 
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L,  but  is  TVD  (or  stable  in  another  relevant  norm)  for  first  order  Euler,  backward  in  time: 

^^"+1  =  (1.10) 

This  can  be  achieved,  for  hyperbolic  conservation  laws,  by  solving  the  backward  in  time 
version  of  (1.2): 

ut  =  f{u),.  (1.11) 

Numerically,  the  only  difference  is  the  change  of  upwind  direction.  Clearly,  L  can  be  com¬ 
puted  with  the  same  cost  as  that  of  computing  L.  We  then  have  the  following  lemma: 


Lemma  1.2.  [12]  The  Runge-Kutta  method  (1.8)  is  TVD  under  the  CFL  coefficient  (1.7): 

Otik 


c  =  mm  ■ 


rrifei’ 

provided  that  a,*  >  0,  and  L  is  replaced  by  L  for  negative  I3ik. 


(1.12) 


Notice  that,  if  for  the  same  k,  both  and  must  be  computed,  the  cost  as 

well  as  storage  requirement  for  this  k  is  doubled.  For  this  reason,  we  would  like  to  avoid 
negative  ^ik  as  much  as  possible.  In  [12],  two  L's  were  used  to  give  a  fourth  order  TVD 
Runge-Kutta  method  with  a  CFL  coefficient  c  =  0.87.  We  will  improve  it  in  this  paper, 
however  unfortunately  we  also  prove  that  all  four  stage,  fourth  order  Runge-Kutta  methods 
with  positive  CFL  coefficient  c  in  (1.12)  must  have  at  least  one  negative 

For  large  scale  scientific  computing  in  three  space  dimensions,  storage  is  usually  a  paramount 
consideration.  Therefore,  there  are  discussions  about  low  storage  Runge-Kutta  methods  [15], 
[1],  which  only  require  2  storage  units  per  ODE  equation.  We  will  consider  in  this  paper 
TVD  properties  among  such  low  storage  Runge-Kutta  methods. 

In  the  next  section,  we  will  give  numerical  evidence  to  show  that,  even  with  a  very  nice 
second  order  TVD  spatial  discretization,  if  the  time  discretization  is  by  a  non-TVD  but 
linearly  stable  Runge-Kutta  method,  the  result  may  be  oscillatory.  Thus  it  would  always  be 
safer  to  use  TVD  Runge-Kutta  methods  for  hyperbolic  problems. 
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The  investigation  of  TVD  time  discretization  can  also  be  carried  out  for  the  general¬ 
ized  Runge-Kutta  methods  (which  have  more  than. one  step)  in,  e.g.,  [6]  and  [7].  We  have 
performed  this  study  but  failed  to  find  good  (in  terms  of  CFL  coefficients  and  whether  L 
appears)  TVD  methods  in  this  class.  The  result  will  not  be  discussed  in  this  paper. 

2  The  Necessity  to  Use  a  TVD  Time  Stepping:  A 
Numerical  Example 


In  this  section  we  will  show  a  numerical  example,  using  the  standard  minmod  based  MUSCL 
second  order  spatial  discretization  [14].  We  will  compare  the  results  of  a  TVD  versus  a 
non-TVD  second  order  Runge-Kutta  time  discretizations.  The  PDE  is  the  simple  Burgers 
equation 

with  a  Riemann  initial  data: 


if  a;  <  0 
if  a:  >  0 


in  (2.1)  is  approximated  by  the  conservative  difference 


(2.2) 


Ax 


where  the  numerical  flux  defined  by 


rr  [fi+i  -  fi-i)  ’ 


fM  =  !> 


with 


u.  1  =  Uj-f— mmmod(Mj+i— Uj, Mj— Uj_i),  —  Uj^i  —  -minmod{uj^2  uj) 


j+j  •'  2 

The  monotone  flux  h  is  the  Godunov  flux 


, ,  _  ,  min„-<„<„+  (4)  ,  if «  < 

h(u  ,u^)  =  <  >„2\  .p  _  ^  4. 

^  (  max„->„>„+  ,  if « 


4 


and  the  now  standard  minmod  function  is  given  by 


minmod{a,  b) 


sign{a)  +  sign{b) 

2 


min(|al,  |6|). 


It  is  easy  to  prove,  by  using  Harten’s  Lemma  [4],  that  the  Euler  forward  time  discretization 
with  this  second  order  MUSCL  spatial  operator  is  TVD  under  the  CFL  condition  (1-6): 


At  < 


Ax 

2  maxj  |u” 


(2.3) 


Thus  At  =  - — will  be  used  in  all  our  calculations. 

^  IXldiXjf  I 

The  TVD  second  order  Runge-Kutta  method  we  consider  is  the  one  given  in  [12]: 


=  w”  +  AfT(u")  (2.4) 

the  non-TVD  method  we  use  is: 

y(i)  =  u”-20AtT(u”)  (2.5) 

=  «"  +  ^AiL(«")  - 

It  is  easy  to  verify  that  both  methods  are  second  order  accurate  in  time. 

If  the  operator  L  is  linear  (for  example  the  first  order  upwind  scheme  applied  to  a  linear 
PDE),  then  both  Runge-Kutta  methods  (actually  all  the  two  stage,  second  order  Runge- 
Kutta  methods)  yield  identical  results  (the  two  stage,  second  order  Runge-Kutta  method  for 
a  linear  ODE  is  unique).  However,  since  our  L  is  nonlinear,  we  may  and  do  observe  different 
results  when  the  two  Runge-Kutta  methods  are  used. 

In  Figure  1  we  show  the  result  of  the  TVD  Runge-Kutta  method  (2.4)  and  the  non-TVD 
method  (2.5),  after  the  shock  moves  about  50  grids  (400  time  steps  for  the  TVD  method, 
528  time  steps  for  the  non-TVD  method).  We  can  clearly  see  that  the  non-TVD  result  is 
oscillatory  (there  is  an  overshoot). 
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Figure  F  Second  order  TVD  MUSCL  spatial  discretization.  Solution  after  500  time  steps. 
Left:  TVD  time  discretization  (2.4);  Right:  non-TVD  time  discretization  (2.5). 


Such  oscillations  are  also  observed  when  the  non-TVD  Runge-Kutta  method  coupled  with 
a  second  order  TVD  MUSCL  spatial  discretization  is  applied  to  a  linear  PDE  (uf  -1-  Ux  0). 
Moreover,  for  some  Runge-Kutta  methods,  if  one  looks  at  the  intermediate  stages,  i.e. 
for  1  ^  i  m  in  (1.8),  one  observes  even  bigger  oscillations.  Such  oscillations  may  render 
difficulties  when  physical  problems  are  solved,  such  as  the  appearance  of  negative  density 
and  pressure  for  Euler  equations  of  gas  dynamics.  On  the  other  hand,  TVD  Runge-Kutta 
method  guarantees  that  each  middle  stage  solution  is  also  TVD. 

This  simple  numerical  test  convinces  us  that  it  is  much  safer  to  use  a  TVD  Runge-Kutta 
method  for  solving  hyperbolic  problems. 

3  The  Optimal  TVD  Runge-Kutta  Methods  of  Sec¬ 
ond,  Third  and  Fourth  Order 

In  this  section  we  will  try  to  identify  the  optimal  (in  the  sense  of  CFL  coefficient  and  the 
cost  incurred  by  L  if  it  appears)  TVD  Runge-Kutta  methods  of  m-stage,  m-th  order,  for 
m  =  2,3,4,  written  in  the  form  (1.8). 

For  second  order  m  =  2,  we  can  choose  /3io  ^nd  Q21  free  parameters.  The  other 


I 


1 
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coefficients  are  then  given  as  [12]: 


'  aio  =  1 
020  =  1  ~  0:21 
/520  =  1  -  ^  -  O2l/?10 

.  ^21  =  2;^ 


(3.1) 


Proposition  3.1.  If  we  require  aik  >  0  and  /?jjt  >  0,  then  the  optimal  second  order  TVD 
Runge-Kutta  method  (1.8)  is  given  by 

ui^)  ^  u^  +  AtL{u^)  (3.2) 


with  a  CFL  coefficient  c  =  1  in  (1.9). 

Proof:  If  we  would  like  a  CFL  coefficient  c>  1,  then  oio  =  1  implies  ^lo  <1,  which  in  turn 
implies  ^  >  I-  Also,  021  >  /^2i  =  which  implies  q;2i/?io  >  |.  We  thus  have 


/?20  —  1  ~ 


1 

2^10 


~  O2i/3i0  ^  ^  ~  2 


which  is  a  contradiction. 


□ 


For  the  third  order  case  m  =  3,  the  general  Runge-Kutta  method  consists  of  a  two 
parameter  family  as  well  as  two  special  cases  of  one  parameter  families  [9] .  We  can  similarly 
prove  the  following  proposition: 

Proposition  3.2.  If  we  require  aik  >  0  and  /3ik  >  0,  then  the  optimal  third  order  TVD 
Runge-Kutta  method  (1.8)  is  given  by 

=  u^  +  AtL{u^)' 

tfC)  =  q.  -t-  ^AtL{u^^^)  (3.3) 

=  iu-  +  ^„(2)  + 

with  a  CFL  coefficient  c  =  1  in  (1.9). 

Proof:  The  proof  is  more  technical,  and  is  given  in  the  Appendix. 
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□ 

For  the  fourth  order  case  m  =  4,  the  general  Runge-Kutta  method  again  consists  of  a  two 
parameter  family  as  well  as  three  special  cases  of  one  parameter  families  [9].  Unfortunately, 
this  time  we  cannot  avoid  the  appearance  of  negative  ^iki 

Proposition  3.3.  The  four  stage,  fourth  order  Runge  Kutta  scheme  (1.8)  with  a  nonzero 
CFL  coefficient  c  in  (1.12)  must  have  at  least  one  negative  Afc. 

Proof:  The  proof  is  technical,  and  is  given  in  the  Appendix. 

□ 


We  thus  must  settle  for  finding  an  efficient  solution  containing  X,  which  maximizes 
where  c  is  the  CFL  coefficient  (1.12)  and  i  is  the  number  of  Ls.  This  way  we  are  looking  for 
a  TVD  method  which  reaches  a  fixed  time  T  with  a  minimal  number  of  residue  evaluations 
for  L  or  L.  We  use  a  computer  program  and  the  help  of  optimization  routines  to  achieve 
this  goal.  The  following  is  the  best  method  we  can  find: 


u”  +  -AtL(n”) 

£ji 


^  _  10890423  ^  ^AtT(u(^)) 

1600  25193600  ^  ^1600  7873  ^ 

5.3989  „  102261  .  ,  4806213  . 

2500000''  5000000^  ^^"  ^^  20000000" 

5121  ^  f  /i\>.  23619  (2)  7873  A,r/ 

g"”  +  +  S""’  +  +  5“"’  + 


with  a  CFL  coefficient  c  =  0.936  in  (1.12).  Notice  that  two  Ls  must  be  computed.  The 
effective  CFL  coefficient,  comparing  with  an  ideal  case  without  Ls,  is  0.936  x  |  =  0.624. 
Since  it  is  difficult  to  solve  the  global  optimization  problem,  we  do  not  claim  that  (3.4)  is 
the  optimal  4  stage,  4th  order  TVD  Runge-Kutta  method. 
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4  The  Low  Storage  TVD  Runge-Kutta  Methods 

For  large  scale  scientific  computing  in  three  space  dimensions,  storage  is  usually  a  paramount 
consideration.  Therefore,  there  are  discussions  about  low  storage  Runge-Kutta  methods  [15], 
[1],  which  only  require  2  storage  units  per  ODE  variable.  We  will  consider  in  this  section 
TVD  properties  among  such  low  storage  Runge-Kutta  methods. 

The  general  low-storage  Runge-Kutta  schemes  can  be  written  in  the  form  [15],  [1]': 

„(«•)  =  +  i  =  l,...,m  (4.1) 

=  u",  Ao  =  0  . 

Only  u  and  du  must  be  stored,  resulting  in  two  storage  units  for  each  variable. 

Carpenter  and  Kennedy  [1]  have  classified  all  the  three  stage,  third  order  (m=3)  low 
storage  Runge-Kutta  methods,  obtaining  the  following  one  parameter  family: 

/36c|  4-  36ci  -  135ci  +  84c2  -  12 
2c2  ~\~  C2  —  2 

124  -  IS4  +  18c^  -  IIC2  -t-  2 
36c2  —  36c2  +  134  —  3c2  -f  4 
69c^  -  62c^  -f  28c2  -  8 
.34c^  -  46c^  +  34c^  -  13c2  -b  2 

C2  (4.2) 

12c2(c2  -  1)(3^2  -  Z-i)  -  (3^2  - 
144C2(3C2  -  2)(c2  -  1)2 
-24(3c2  -  2)(c2  -  1)" 

(.32:2  -  ^1)^  -  12c2(c2  -  1)(322  -  ^l) 

-Zi(6C2  -  4C2  +  1)  +  3z3 
(2c2  -b  1)21  -  3(c2  +  2)(2c2  -  1)2 
— Z4Z1  -f-  108(2c2  —  1)^2  —  3(2c2  —  1)25 

242iC2(c2  —  1)"*  -|-  72C226  “b  72c2(2c2  —  13) 


Zi  = 

z-i  = 
Z3  = 

Z4  = 

Z5  = 
26  = 
By  = 
B2  = 

Bs  = 
A2  — 
As  = 
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Which  is  converted  into  the  form  of  equation  (1.8),  by  introducing  three  new  parameters. 
Then  we  search  for  values  of  these  parameters  that  would  maximize  the  CFL  restriction, 
again  by  a  computer  program.  The  result  seems  to  indicate  that 

C2  =  0.924574  (4.3) 

gives  an  almost  best  choice,  with  CFL  coefficient  c  =  0.32  in  (1.9).  This  is  of  course  less 
optimal  than  (3.3)  in  terms  of  CFL  coefficients,  however  the  low  storage  form  is  useful  for 
large  scale  calculations. 

Carpenter  and  Kennedy  [1]  have  also  given  classes  of  5  stage,  4th  order  low  storage  Runge- 
Kutta  methods.  We  have  been  unable  to  find  TVD  methods  in  that  class  with  positive  aik 
and  ^ik.  Notice  that  L  cannot  be  used  without  destroying  the  low  storage  property,  hence 
negative  l3ik  cannot  be  used  here. 

5  Concluding  Remarks 

We  have  given  a  simple  but  illustrating  numerical  example  to  show  that  it  is  in  general  much 
safer  to  use  a  TVD  Runge-Kutta  method  for  hyperbolic  problems.  We  then  explore  the 
optimal  second,  third  and  fourth  order  TVD  Runge-Kutta  methods.  While  for  second  and 
third  order  optimal  methods  are  found  with  a  CFL  coefficient  equal  to  one,  for  fourth  order 
we  simply  give  the  best  method  we  can  find.  A  TVD  third  order  low  storage  Runge-Kutta 
method  is  found,  which  uses  only  two  storage  units  per  equation  and  has  a  CFL  coefficient 
equal  to  0.32.  Finally,  we  prove  that  general  four  stage  fourth  order  Runge-Kutta  methods 
can  not  be  TVD  without  introducing  an  adjoint  operator  L. 

Acknowledgments:  We  would  like  to  thank  Mordechai  Berger  and  Mark  Carpenter  for 
helpful  discussions. 
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6  Appendix 

In  this  appendix  we  prove  Proposition  3.2  and  Proposition  3.3. 

We  write  the  general  4  stage,  4th  order  Runge-Kutta  method  in  the  following  standard 

form  [9]: 

=  u”  -f  cioZ/(u”) 

=  u”  -|-  C2oAtiy(u”)  -|-  C21 

t<(3)  =  U"  +  c3oAtL(u”)  +  C31  +  C32AtL(u(''))  (6.1) 

=  u”  -f  C4oAtX(u")  -1-  C4iAtL(«^^^)  -I-  C42AtL{u^'^^)  ->I-  C4zAtL{u^‘^^) 

The  relationship  between  the  coefficients  Cik  here  and  aik  and  ffk  in  (1.8)  is: 

ClO  =  /^lO 

C20  =  ^20  +  0^21  ^10 

C21  =  ^21 

C30  =  0:320:21  +  Q!31^10  +  «32/?20  + /^SO 


12 


(6.2) 


C31  =  0!32^21  +  /?31 

C32  =  032 

C40  =  Q;430!320:21/^10  +  Q!43Q!32^20  +  Q;43Q!31^10  +  Q!420!21^10 

+0:41^10  +  Oii202O  +  OCA303O  +  0iO 

C41  =  a43<^32021  +  0142021  +  0-43031  +  041 

C42  =  043032  +  042 

C43  =  043 

For  a  third  order  Runge-Kutta  method,  the  general  form  (6.1)  is  similar  without  the  last 
line  (and  with  replaced  by  The  relationship  (6.2)  also  is  similar  without  the  last 

four  lines  for  C40,  C41,  C42  and  C43. 

Proof  of  Proposition  3.2: 

The  general  third  order,  three  stage  Runge-Kutta  method  in  the  form  (6.1)  is  given  by  a 
two  parameter  family  as  well  as  by  two  special  cases  of  one  parameter  families  [9]. 

•  General  Case:  If  03  ^  02,  03  /  0,  02  7^  0  and  02  7^  |: 

=  02 

_  .80203(1  -  02)  -  ol 
02(2  -  .802) 

_  03(03  -  02) 

02(2  —  802) 

_  j  2  —  8(02  -\-  03) 

602O3 

_  803  —  2 

602(03  —  02) 

2  —  802 
603(03  -  02) 

Notice  that  602C21C32  =  1  and  C20  +  C21  =  03.  If  we  want  to  have  a  CFL  coefficient 
c  >  1  in  (1.9),  we  would  need  0,^  >  Afc  >  0  unless  both  of  them  are  zeroes.  This  also 


Cio 

C20 

C21 

C30 

C31 

C32 
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implies  that  dk  >  0  by  (6.2).  Also,  notice  that  =  A,i-i  >  0,  otherwise  that  stage 
is  not  necessary. 

Now,  Cio  =  Ao  <  <^10  =  1  ^10  >  0  imply  0  <  0!2  <  1. 


1.  as  >0:2. 

C21  >  0  implies  02  <  |,  and  C31  >  0  requires  03  >  |. 

^20  >  0  and  021  >  /^2i  imply  C20  >  O2i/?io  >  /?2iAo,  which  is  03  -  C21  >  02102,  or 
>  C21.  So  we  must  have 

1+02 

302  —  2o2 

03  <  — — - 

1  +  0:2 

On  the  other  hand,  1S31  >  0  requires  C31  >  032^21  >  C32C21  =  which  is 
3o3  -  2  >  03  -  0:2,  or 

0:3  >  1  -  2*^2 

Combining  these  two  inequalities,  we  get  1  —  |a2  <  or  (2  -  3o2)(l  -  0:2)  < 

0,  which  is  a  contradiction,  since  2  —  3o2  >  0  and  1  —  02  >  0. 


2.  02  >  0:3. 


0:3  =  C20  +  C21  >  0  requires  03  >  0. 

C32  >  0  requires  0:2  >  |,  and  C31  >  0  requires  0:3  <  |. 
C31  >  ^32/321  >  C32C21  =  6^,  which  is 


03  <  1  -  2^2, 


C20  >  021^10  >  /^2n^io  requires 


03  > 


02(3  —  202) 
1  +  02 


Putting  these  two  inequalities  together,  we  have  <  1  —  2^25  which  means 

(2  —  3o2)(1  —  0:2)  >  0,  a  contradiction  since  1  —  02  >  0  and  2  —  3o2  <  0. 
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•  Special  Case  I:  02  =  eta  =  §•  In  this  case 

2 
3 

2  1 

3  AtJs 
1 

1 

4 

3 

—  —  U}j 

4 

C32  =  U>3 

i^3i  ^  0  &nd  a32  >  ^32  =  C32  requires  C31  >  0132/^21  >  Cz2^2\  =  4  which  implies  tua  <  |. 

^20  >  0  and  021  >  ^21  =  C21  requires  C20  >  0:21^10  >  |c2i,  which  means  |  -  >  §5^, 

for  which  we  must  have  wa  >  |.  A  contradiction. 

•  Special  Case  II:  03  =  0.  In  this  case  the  equations  read 

2 
3 

1 

dwa 
_ 1_ 

4a;3 

1 

--a;3 

3 

4 

^3 

Clearly  C20  and  C21  cannot  be  simultaneously  nonnegative. 

•  Special  Case  III:  02  =  0.  In  this  case  the  method  is  not  third  order. 

■  □ 


Cio  — 
C20  = 

C21  = 

C30  = 
C31  = 
C32  = 


ClO  — 

C20  = 

C21  = 

C30  = 
C31  = 
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Proof  of  Proposition  3.3: 

Recall  that  all  the  aik's  must  be  nonnegative  to  satisfy  our  TVD  criteria.  From  the 
relationship  (6.2)  between  the  coefficients  of  (6.1)  and  of  (1.8),  we  can  see  that  nonnegative 
Pik's  imply  nonnegative  Ci^s.  We  now  show  that  we  cannot  have  all  nonnegative  cafe’s. 

•  General  Case.  If  two  parameters  02  and  03  are  such  that:  02  ^  03,  02  /  1,  02  0, 

^2  /  tts  /  1,  «3  ^  0,  03  ^  5,  and  602^3  -  4(02  +  0C3)  +  3^0.  Then  the  coefficients 
dk  are  [9]: 

_  03(03—02) 

Cio  =  012,  C20  —  O3  —  C21,  C21  —  202(1-202)' 

.  _  (l-a2)[a2+a3-l-(2a3-lffi  ^  _  (l-2a2)(l-a2)(l-a3) 

C30  —  i  C31  C32,  C31  —  2a2(a3-a2)l6a2^<3-4(or2+o;3)+3]  ’  Q3(a'3— o:2)[6<>20;3-4(a2+«3)+3]  ’ 

1  ,  1-2(024-03)  _ _ 2oa-l  „  _ _ (1-202)  .  _  1  I  2(o2-|-03)-3 

C40  2  120203  5^41  1202(o3-02)(1— 02)  ’  1203(03-02  )(l-03)  ’  2  12(1— 02)(1— 03) 

There  are  five  possibilities  to  consider: 

1.  02  <  0  implies  c\q  <  0. 

2.  03  >  02  >  0  and  0  <  02  <  |: 

C41  >  0  requires  03  >  |.  C20  >  0  requires  03  <  802  —  4o2  ^  C32  >  0  and 

C31  >  0  require  that  02  >  2  -  603  +  4o|.  Since  this  is  a  decreasing  function  of  03 
when  03  <  ^,  we  obtain  02  >  2  -  5(3o2  -  4o|)  +  4(3o2  -  ialf.  Rearranging, 
we  find  that  0  >  2((2o2  -  1)^  +  4o^)  (2o2  -  1)^,  which  is  impossible. 

3.  03  <  02  and  02  >  |: 

C42  >  0  requires  0  <  03  <  1. 

We  can  only  have  C32  >  0  in  one  of  two  ways: 

(a)  If  1  -  02  >  0,  and  602O3  -  4(o2  +  03)  +  3  >  0. 

C41  >  0  requires  03  <  |.  Simple  calculation  yields  C30  =  1  -  C31  -  C32  = 

(2-6o2+4o^)4-(-5-^i5o2-i2oi)o3-f(4-i2o2ti2o|)o|  ^  q  requires 

2o'2Q;3(6o20f3""4(a2+0(3)+3) 

A+Bas+Cal  =  (2-6o2+4o^)+(-54-15o2-12o^)o3+(4-12o2+12o^)o|  >  0 
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It  is  easy  to  show  that,  when  |  <  a2  <  1,  we  have  ^4  <  0,  <  0  and  (7  >  0. 

Thus,  for  0  <  as  <  |,  we  have 

A+Baz+Cal  <  max  A  + 

which  is  a  contradiction. 

(b)  02  >  1,  and  602O3  —  4(o;2  +  as)  +  3  <  0. 

C31  >  0  requires  Q2  +  os  —  1  —  (2a3  —  1)^  <  0,  which  implies 

(1  -  4a3)(l  -  ctz)  =  4a2  _  5^3  +  1  >  ^2  -  1  >  0. 

Clearly,  this  is  true  only  if  0:3  < 

Now,  C40  >  0  requires  that  0  <  —  2{a2  A  0:3)  +  1  =  2az{Za2  —  1)  +  (1  — 

202)  <  |(3a2  ~  1)  +  (1  “  202)  =  |(1  —  ^2),  an  apparent  contradiction. 

4.  0  <  a2  <  I  and  az  <  ct2'  in  this  case  we  can  see  immediately  that  C42  <  0. 

5.  If  I  <  02  <  0:3,  C21  <  0. 

•  If  6o2Q!3  —  4(0:2  +  otz)  +  3  =  0,  or  if  02  =  0  or  if  03  =  1,  then  this  method  is  not  fourth 
order  [9]. 

•  Special  Case  I.  If  02  =  03  the  method  can  be  fourth  order  only  if  02  =  03  =  |.  In 

this  case  [9]  cio  =  C20  =  |  -  C21  =  C30  =  0,  C31  =  1  -  3wz,  C32  =  dwz, 

C40  —  g,  C41  —  2  ^3^  *^42  ^^3>  ^43  0- 

Clearly  we  need  to  have  C42  =  103  >  0.  To  have  C31  =  1  —  3u;3  >  0  and  C20  =  |  ^  >  0, 

we  require  1^3  =  |.  This  leads  to  the  classical  fourth  order  Runge-Kutta  method. 
Clearly,  then,  021  =  =  —2^20-  This  is  only  acceptable  if  021  =  /?2o  =  0.  But 

021  =  so  in  the  case  where  all  0ikS  are  nonnegative,  the  CFL  coefficient  (1.12)  is 
equal  to  zero. 

•  Special  Case  II.  If  02  =  1,  the  method  can  be  fourth  order  only  if  0:3  =  |.  Then 
[9]  cio  =  1,  C20  =  |,  C21  =  |,  C30  =  1  —  C31  —  C32,  C31  =  — C32  = 

C44  —  g  ^4?  ^42  3?  ^43  ^4* 


+  ic)  =  max  (a,  i(l  -  2a2)(l  -  ^2))  <  0, 


1  „  _  1 

3wi'>  ^“*0  —  0) 


In  this  case  we  want  C31  =  -3^  >  0  which  means  W4  <  0.  But  then  C43  -  u;4  <  0.  So 
this  case  does  not  allow  all  nonnegative  A^’s. 


•  Special  Case  III.  If  03  =  0  the  method  can  be  fourth  order  only  if  0:2  =  Then  [9] 
Cio  =  5,  C20  =  C21  =  C30  =  1  -  C3I  -  C32,  C31  =  |,  C32  =  6w3,  C40  =  I  - 


C4I  =  h  C42  =  'W^3,  C43  — 


Clearly,  C20  =  =  — C21,  one  of  these  must  be  negative.  Thus,  this  case  does  not 


allow  all  nonnegative  life’s,  either. 


□ 
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